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Abstract 

The recent discovery of active brown fat in liuman adults lias led to renewed interest in the role of this key metabolic tissue. 
This is particularly true for neurodegenerative conditions like Huntington disease (HD), an adult-onset heritable disorder 
with a prominent energy deficit phenotype. Current methods for imaging brown adipose tissue (BAT) are in limited use 
because they are equipment-wise demanding and often prohibitively expensive. This prompted us to explore how a 
standard MRI set-up can be modified to visualize BAT in situ by taking advantage of its characteristic fat/water content ratio 
to differentiate it from surrounding white fat. We present a modified IVIRI protocol for use on an 11.7 T small animal MRI 
scanner to visualize and quantify BAT in wild-type and disease model laboratory mice. In this application study using the 
R6/2 transgenic mouse model of HD we demonstrate a significantly reduced BAT volume in HD mice vs. matched controls 
(n = 5 per group). This finding provides a plausible structural explanation for the previously described temperature 
phenotype of HD mice and underscores the significance of peripheral tissue pathology for the HD phenotype. On a more 
general level, the results demonstrate the feasibility of MR-based BAT imaging in rodents and open the path towards 
transferring this imaging approach to human patients. Future studies are needed to determine if this method can be used 
to track disease progression in HD and other disease entities associated with BAT abnormalities, including metabolic 
conditions such as obesity, cachexia, and diabetes. 
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Introduction 

Brown adipose tissue (BAT) is a key tissue for the regulation of 
wliole-body energy metabolism [I]. Ttirougti ttie expression of 
mitochondrial uncoupling protein I (UCP-I), which allows proton 
leakage across the inner mitochondrial membrane, BAT is capable 
of dissipating energy in the form of heat. BAT is thus the principal 
effector tissue of adaptive, non-shivering thermogenesis. A series of 
breakthrough positron emission tomography (PET) and computed 
tomography (CT) imaging studies has recently shown BAT 
deposits even in adult humans [2-5]. This resulted in a renewed 
interest in the role of BAT especially for metabolic conditions such 
as obesity and cachexia, thereby creating the need for appropriate 
imaging technologies [6]. Large-scale studies aiming at the 
structural and functional characterization of BAT by PET/CT 
in human adults are still scarce [4,6,7]. Factors contributing to this 
limitation include the high costs for the required tracers, safety 



concerns about the radiation involved, and the required activation 
of BAT either via cold-exposure or pharmacologically. 

Huntington disease (HD), an autosomal dominant adult-onset 
neurodegenerative disease caused by the abnormal expansion of 
the CAG-repeat tract in the huntingtin gene, represents an area of 
research where these new developments are of particular interest 
[8] . In addition to the classic triad of neuro-psychiatric symptoms, 
i.e. movement disorder, cognitive decline and psychiatric abnor- 
malities, HD patients display a progressive energy deficit that is of 
high clinical and scientific relevance [9]. This energy deficit is 
recapitulated in transgenic mouse models of HD and this led to the 
implication of brown fat dysfunction in HD pathogenesis, even 
before BAT-activity in human adults was discovered [10,11]. 
Selected transgenic mouse models of HD therefore provide a 
useful model system to study disease-related BAT dysfunction 
in vivo. 

Hamilton and colleagues performed an in-depth analysis of the 
magnetic resonance (MR) signatures of BAT and white adipose 
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tissue (WAT) [12]. They reported several key MR properties tiiat 
are different between BAT and WAT, including fat fraction, 
changes in the Tl relaxation rate of the water-bond protons, and 
the degree of lipid saturation. The use of the histological and 
physiological properties of BAT has the major advantage of not 
requiring prior stimulation and has been applied for fat 
characterization by means of magnetic resonance spectroscopy 
(MRS) and chemical-shift imaging (CSI) [13-15]. 

With the exception of intermolecular zero-quantum coherences 
spectroscopy aiming for assessment of the intravoxel distribution of 
fat and water spins [16] and blood oxygen level dependent MRI 
[17], the majority of published fat fraction imaging methods by 
MRI relies on direct fat/water quantification as initially 
introduced by Dixon and colleagues in 1984 [18]. Based on the 
original publication, several modifications of the technique have 
been proposed [19]. Although in principle off-resonance contri- 
butions can be compensated by the multi-point techniques, 
additional non-static phase errors and the need for phase 
unwrapping quite frequently introduce errors in the resulting 
fat/water maps and may cause exclusion of data [15] especially at 
ultrahigh field strength. 

The goal of this work was to establish and validate a simplified 
two-point technique for direct visualization of BAT, using the 
difference between the magnitude of images obtained with in- 
phase and opposed-phase condition. The applicability of the 
proposed technique was tested pre-clinically in an HD disease 
model (R6/2 transgenic mouse). HD appears of special interest 
since brown fat abnormafities are of emerging interest in the 
investigation of the HD energy deficit, and the clear genetic 
pathogenesis of HD provides an excellent basis for future 
mechanistic studies in transgenic model systems and facilitates 
validation in human patients and translational approaches. 

Results 

All scans could be completed successfully in the investigated 
cohorts. Automatic thresholding applying the Otsu method was 
effective in full suppression of background signal as well as pixels 
with low signal intensities as e.g. present in the lung and cortical 
bone. Figure 1 shows a direct in vivo comparison of the 2PM 
approach with the conventional 3PD method. The image intensity 
in the 2PM approach is linearly related to the fat-water fraction 
with a maximum at equal magnetization of fat- and water spins. 
Since BAT is expected to show an almost equal magnetization of 
fat- and water spins, BAT presents as bright signal intensity in the 
final 2PM images. Water- and fat-only voxels show a clearly 
reduced intensity. The signal in the fat fraction images derived 
from the 3PD method scales linearly with the fat fraction in the 
magnetization. Water-only voxels appear dark, fat-only voxels 
appear bright, and BAT voxels present with intermediate signal 
intensity. The direct comparison shows an excellent resemblance 
of the resulting BAT distribution between the two approaches. 
The interscapular fat deposits (Figure 1, arrows) can be clearly 
appreciated as well as some smaller pads of BAT extending 
anteriorly and bilaterally around a large vessel (Figure 1, 
arrowheads). Note that both methods show a similar sensitivity 
to partial volume effects. 

An overview of the BAT pattern for the wild-type and 
transgenic R6/2 mice is provided in Figure 2 for a slice centered 
at the location of tlu' prominent interscapular BAT. As expected 
and as described previously [20,21], the reduced body size and 
general organ atrophy of R6/2 mice vs. wild-type controls is 
immediately evident. 



Even though the spatial resolution was lowered for reduction of 
the acquisition time, the BAT distribution can be readily retrieved 
from the resulting 2PM images. Semi-automatic quantification of 
the BAT volumes showed a highly significant positive correlation 
of the 2PM and the conventional 3PD method (R = 0.976, 
p<0.01, Spearmans's rho). Further analysis allowed for direct and 
quantitative comparison of BAT volumes in the region of interest 
between HD mice and matched wUd-type controls. As shown 
in Figure 3a, R6/2 mice have a greatly reduced BAT volume 
(P<0.01) of 18.46 nl (min: 14. 18 [ll, max: 23.85 ml vs. 36.82 1x1 
(min: 32.34 pil, max 45.40 |il) in controls. The difference is blunted 
but stiU manifest and statistically significant (P<0.05) when 
normalization to total body volume as a correction for the 
different body sizes is introduced (Figure 3b). 

Discussion 

The study e\'aluates a non-invasive method, which allows for 
in vivo and ex vivo BAT volume quantification from magnitude 
images. The intrinsic insensitivity of the technique to any kind of 
phase errors introduced by system imperfections should cause an 
at least similar reproducibility of 2PM as the more established fat- 
fraction methods [14,22,23]. Since the approach is based on a 
conventional FLASH technique, it can be readily used on any 
conventional MRI scanner without further modification of 
acquisition methods. In comparison to other MRI-based BAT- 
imaging techniques, the advantage of the 2PM method is the 
straightforward approach, which does not require any special 
expertise in advanced MRI techniques or image processing. The 
technirjue can be readily modified further to provide blood oxygen 
level dependent contrast. Khanna et al. [17] showed that T2* 
weighting times of about 2 ms are sufficient for the assessment of 
BAT activity. For magnetic field strengths between 1.5 and 1 1.7 T 
as used for clinical and pre-clinical imaging, one of the required in- 
and opposed-phase echo times can be chosen close to the 
suggested value. Furthermore, the 2PM method can also easily 
be adapted to the suggested asymmetric spin echo acquisition if 
the required echo times cannot be met due to gradient power 
constraints. In this context the possibility of intrinsic identification 
of the BAT location from the scans obtained at normal body 
temperature is advantageous. Compared to the more complex 
intermolecular zero-quantum MRI [16,24], 2PM is hampered by 
partial volume effects similar to the fat fraction methods. However, 
slightiy higher spatial resolution can be realized without significant 
scan time penalties since images with only two different echo times 
are required for providing the BAT information. 

By use of our MRI-based tissue analysis, we show that R6/2 
mice have significantly reduced BAT, even when corrected for 
reduced overall body size. This finding provides a plausible 
structural explanation for the striking cold-sensitive phenotype of 
R6/2 mice and other HD models [10,1 1]. The results are also in 
accordance with earlier reports of structural BAT abnormalities in 
transgenic mouse models of HD, including the R6/2 strain, all of 
which were obtained through non-survival conventional tissue 
dissection [10,11,25,26]. Importantiy, our MRI-method now 
provides a practical strategy for monitoring BAT structure in vivo 
and longitudinally in the same individual animal. 

While we focused our present study on the R6/2 mice, other 
transgenic mouse strains with a cold-sensitive body temperature 
phenotype such as the N171-82Q, model of HD [10] and the 
mutant SODl model of amyotrophic lateral sclerosis [27] should 
also be investigated with respect to BAT structure, distribution and 
function in vivo. We would like to emphasize that this approach 
can naturally be applied to other diseases, e.g. obesity and 
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Figure 1. Direct comparison of the novel two-point magnitude method (2PM) and the conventional three-point phase-based Dixon 
method (3PD) in a wild-type mouse. Interscapular BAT deposits (arrows) as well as smaller pads of BAT (arrowheads) can be appreciated. Please 
note the excellent agreement between the two methods. Due to the suppression of pure fat and water voxels in the 2PM technique, a clear 
improvement of contrast between WAT and BAT areas can be appreciated with the 2PiVl technique. 
doi:1 0.1 371 /journal.pone.01 05556.g001 



An important limitation of the data set presented here is the 
focus on a single fat pad in the interscapular region. We selected 
this region as it is recognized as the largest and most prominent 
BAT depot in rodents [28]. This choice permitted reliable 
anatomical localisation of the region of interest. It also allowed 
for direct comparison with the histological data from HD models, 
as to the best of our knowledge all the BAT tissue dissection studies 
in HD used interscapular BAT. This focus, however, did preclude 
characterisation of the BAT distribution patterns. As we know 
from imaging studies, fat distribution can vary in specific disease 
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diabetes. The results further suggest that this non-invasive 
methodology can be readily applied to human patients and 
control subjects. However, the dilferent distribution of the BAT in 
humans especially in small pads as e.g. as in perivascular tissue 
may in general limit the direct BAT imaging approaches 
independent of the imaging technique used. In this context the 
improved scan efficiency of the proposed 2PM method may be 
beneficial due to better spatial resolution without sacrificing scan 
time compared to multi-point Dixon methods. 



Figure 2. Comparison of morphological images and resulting BAT-weighed images from R6/2-transgenic (TG) and corresponding 
wild-type (WT) mice. Note the overall reduced tissue volume in the R6/2 mice. The slices were selected to comprise the shoulder region where the 
prominent interscapular BAT pad is located. 
doi:1 0.1 371/journal.pone.01 05556.g002 
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Figure 3. Total BAT volumes a) and BAT volumes normalized to total body volume b) for female HD-transgenic (R6/2) and matched 
wild-type control mice. (n = 5 per group; *P<0.05, **P<0.01; unrelated Wilcoxon Ranl< sum test; columns (error bars) represent means (standard 
deviation)). 

doi:1 0.1 371 /journal.pone.01 05556.g003 
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Figure 4. Resulting magnetization M (black) in a voxel containing almost equal concentration of water- (M„, blue) and fat- (Mf, red) 
bond protons for different phase differences. Respective iVlR images of the mouse neck in axial orientation are exemplarily shown for in-phase 
(left) and opposed-phase (right) condition. The reduction of signal intensity in the interscapular BAT (arrows) is clearly appreciable. 
doi:1 0.1 371/journal.pone.01 05556.g004 
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Figure 5. Comparison of the fat-fraction method based on a three-point Dixon approach (3PD, a), and the proposed two-point 
magnitude based method (2PIVI, b). The fat-fraction a) and BAT b) image are normalized from 0 to 1 . The grey bars indicate the respective grey 
values used as threshold for WAT (>0.8, 3PD) and for BAT (>0.4 and <0.8, 3PD; >0.7, 2PM). WAT cannot be identified by the 2PM method, since 
signal from voxels containing only fat or water will be cancelled out during the involved subtraction. 
doi:10.1371/journal.pone.0105556.g005 



states [29] and how such distribution changes play out in BAT 
pathology remains to be investigated. As a further limitation, our 
study was not designed to detect the so called beige or brite 
(irown-in-whzte) adipocytes [30] but in principle the presented 
tissue characterisation algorithm can be utilized to do so, if 
sufficient spatial resolution was achieved. Further in vivo studies 
with cold exposure have to be performed since the assessment of 
activated BAT has not been tested in this study. Also, since the 
method presented here relies on the characteristic relative fat/ 
water content of BAT the use under conditions where the BAT 
composition is altered, e.g. starvation, require adjustments that 
were beyond the scope of this study. 

In summary, we present a straightforward technique for 
assessing BAT structure with conventional MRI set-ups. Further 
research and development is necessary to integrate functional 
information, e.g. in the form of blood oxygen level dependent 
signal, into the MRI data acquisition protocol. The non-invasive 
in vivo character of the methodology allows longitudinal studies 
that were not previously possible. In a first application of the new 
methodology, we were able to show that symptomatic R6/2 mice 
have a reduced BAT volume. This finding provides insight into to 
the energy deficit in HD and highlights the importance of the 
peripheral tissue pathology for the understanding of the HD 
phenotype [31]. 

Materials and Methods 

Theory 

Depending on the environment, diflFerendy bond 'H protons 
experience different resonance frequencies caused by the shielding 
effects of the electrons. This characteristic chemical shift can be 
utilized for tissue characterization and structural analysis by MRS 



and is frequently applied to chemical analysis [32-34] in 
biomedical research. For H proton spectroscopy, the reference 
with 0 ppm chemical shift is defined by tetramethysilane (TMS) 
and the chemical shift of different organic compounds is defined 
relative to TMS. Chemical shifts of fat bond H protons range 
from 0.9 ppm (-(CHajn-CHj) to 5.3 ppm (-CH = CH-) and 
water (H2O) showing a chemical shift of 4.7 ppm. 

Hamilton et al. [12] investigated the MR properties of WAT 
and BAT. Important findings included alteration of endogenous 
biochemical characteristics as water - fat fraction (wff = ^"T'' , 

with Snriier being the MR signal from water-bond protons and Sfat 
the respective MR signal from fat-bond protons), Tl relaxation 
times, and lipid saturation. From the investigated parameters wff 
showed the most pronounced difference between WAT 
{wff < 0.2) and BAT (wff-l). The higher amount of water in 
BAT has driven the application of fat fraction quantification to 
BAT quantification in mice [14] as well as in humans [15,22,35]. 

Fat fraction quantification demands accurate knowledge of the 
fat and water content in each voxel. Separation of water and fat 
can be achieved by Dixon techniques [18,23,36]. Dixon methods 
rely on the phase difiFerences between MRI images obtained with 
different echo times (TE). The phase differences are introduced by 
the chemical shift of the different compounds. The MRI signal 
5'(?)from voxels containing different compounds results as 

with Mj ( t) being the magnetization of the J-th component at time t 
and A ff'-^ it's chemical shift. Although, different proton moieties 
are present, the fat contribution in BAT and WAT is dominated 
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by the methylene protons (1.3 ppm), followed by methyl (0.9 ppm) 
and allylic protons (2.0 ppm). Considering water and methylene 
protons, only, Eq. 1 simplifies to 

S{t)=M„{t)+Mf{t)e"^('), (2) 

with (p{t) = 2n^fl~^ t being the phase between water and fat-bond 
protons. Theoretically, the fat {Mf) and water {M„) image can be 
easily reconstructed from the sum {My,) and the difference {Mf) of 
gradient echo images acquired in in-phase (p = 2nn) and opposed- 
phase (p = (2m— 1)71) condition as suggested as the original two- 
point Dixon method [18]. A major limitation for the two-point 
Dbcon method arises from local off-resonances, introducing an 
additional TE-dependent phase term and Eq. 2 must be rewritten 
as: 

S{t) = {M^{t) + Mf{t)e'^<''^)e'''v'-'\ (3) 

with A«?(?) = 27tA/'^^? being the off- resonance introduced phase 
error at time t. To compensate for the additional phase term, off- 
resonance maps and phase unwrapping techniques have to be 
used. Several techniques have been introduced so far for 
compensation of the off-resonance induced phase errors [36-38]. 
AH techniques require additional measurements, which prolong 
the overall acquisition time. Off-resonance contributions can be 
compensated by the multi-point techniques, additional non-static 
phase errors and the need for phase unwrapping frequently 
introduce errors in the resulting fat/ water maps especially at 
ultrahigh field strength. 

In this study, a magnitude-based approach is introduced 
suggesting the use of the difference of the magnitude of two 
images acquired with echo times (TE) values for BAT volume 
quantification. The approach is based on the assumption that a 
major difference between WAT, normal biological tissue, and 
BAT results from the almost equal concentration of water- and fat- 
bond protons in BAT. The magnitude of 5'(;) results from [Eq. 3] 
as 

\S{t)\ = \M,,{t) + Mf{t)e'i''^'^, (4) 

where the phase (p (?) between the water- and fat-bond protons 
cause an intensity modulation depending on the magnetization of 
the fat Mf and water M,^ with almost complete signal cancellation 
at opposed phase condition in case of equal water and fat 
contributions as shown in Figure 4. For background suppression 
of voxels containing a single tissue, only, the final BAT signal is 
calculated according to: 

Seat — \Sin—phase\ {^opposed— phasel (5) 

Animal groups 

The suggested two-point magnitude (2PM) technique was 
initially evaluated in vivo in direct comparison with the three- 
point Dixon (3PD) technique in a group of three wild-type mice 
(C57/B6, 12 months, female). All in vivo data acquisition was 
carried out under isoflurane anesthesia (3'X) for induction and 
1.5% for maintenance). 

ApplicabiUty to cohort studies was evaluated in the R6/2 
transgenic mouse model compared to sex-matched wild-type litter 
mates. The R6/2 mouse is among the most commonly used 
transgenic models of HD [20,39]. They display cold-intolerance 
highly suggestive of BAT dysfunction, and structural abnormalities 



of white and brown adipose tissue are well-documented 
[10,25,40]. Female mice were used at 12 weeks, an age when 
this model is clearly symptomatic. Assuming a volume change in 
the order of 50% in the R6/2 group, to achieve a statistical power 
of 80% sample sizes of 5 animals per group were required. Power 
was calculated with nQuery (version 10). To exclude any impact of 
respiratory and cardiac motion artifacts in this feasibility study, the 
animals were sacrificed (CO2 inhalation) immediately prior to the 
MRI investigation. 

In vivo animal experiments were performed in accordance with 
German animal protection laws and had been approved by the 
national animal board (TVA 1001, "Etablierung KleintierbUdge- 
bung an 3T und 11. 7T", Regierungsprasidium Baden- Wurttem- 
berg, Ttibingen, Germany). 

MRI protocol 

The proposed two-point magnitude method (2PM) was 
implemented and evaluated against a 3PD method based on data 
acquired on a dedicated high-field 1 1 . 7T small animal MRI 
system (BioSpec 117/16, Bruker, Germany, Ettiingen). All data 
were acquired applying a 4cm quadrature volume resonator. 

For initial in vivo evaluation (C57/B6) of the 2PM methods, a 
high-spatial resolution spoiled fast low-angle shot gradient echo 
technique (HR-FLASH, Ar= 100x100x500 |xm'^) was applied. 
Data were acquired at five different echo times (TE = 1.203 ms, 
1.313 ms, 1.459 ms, 1.507 ms, 1.751 ms), yielding phase shifts of 
-71/6, 7i/2, 7 71/6, 71, 2 71 respectively. The applicability to cohort 
studies (R6/2, wild-type) was tested with a reduced spatial 
resolution protocol (LR-FLASH, Ar= 120x200x500 tmi^) to 
ensure rapid quantification of the BAT volumes. In all animals, 
the complete neck area was c()\c-rcd by acquiring 9 slices of 
500 |xm slice thickness with 1 mm spacing. Acr|uisition time Tacq, 
for the acquisition of a single TE resulted to T \(;q^~ 3m: 50s (LR- 
FLASH) and Tagq.= 6m 24s (HR-FLASH). Excitation angle a 
and repetition time TR were as a = 15° and TR= 75 ms. 

Data processing 

After acquisition, the data was transferred to an independent 
workstation and further processed for BAT quantification. Fat 
fraction images were derived from the images acquired with -7t/6, 
71/2, 771/6 phase shifts applying a multi-point Dixon method. For 
ensuring high qualit)' of the fat fraction maps, the "fat-water 
toolbox" (http://ismrm.org/workshops/FatWaterl2/ data.htm) 
provided by the ISMRM was apphed for all calculations. For 
evaluation of the suggested two-point magnitude technique, the 
difference of the magnitude of the images obtained with in-phase 
and opposed-phase condition was calculated. An overview of the 
two approaches is provided in Figure 5. For noise reduction prior 
to further analysis, a 3x3 point medium filter was applied. 
Background signal was removed by thresholding the input data. 
The respective threshold value was derived from the in-phase 
images according to the Otsu threshold algorithm [41]. 

Volume quantification of the BAT was performed using a semi- 
automatic segmentation tool (http://www.itksnap.org/, [42]). Re- 
gions of BAT were identified in a slice-by-slice fashion according 
to the mean normalized grey value / with />0.7 for the 2PM 
method and 0.4<7<0.8 for the 3PD method. Falsely identified 
regions due to partial volume effects were removed manually for 
both approaches. The evaluation was performed by an experi- 
enced MRI scientist (VR, > 15 yrs in research). Quantification was 
performed blinded for the two investigated approaches. 

After semi-automated classification, resulting BAT volumes 
were compared between 2PM and 3PD as well as between the 
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HD-transgenic mice and the wild-type controls. Correlation of the 
volumes derived from the investigated methods was assessed 
applying a Spearman's rho test and the significance of BAT 
volume changes between the different groups was assessed by an 
unrelated Wilcoxon rank sum test. P-values below 0.05 were 
considered significant. All statistical analysis was performed with 
SPSS 21.0 (IBM Corp, Armonk, NY). 
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